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We measure the time-dependent CP asymmetry in B^^ — > ii's7r'^7 decays for two regions of K^-tv" 
invariant mass, m{Ks-iv^), using the final BABAR data set of 467 x 10^ BB pairs collected at the 
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—0.36 ± 0.33 ± 0.04. The uncertainties are statistical and systematic, respectively. 
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The radiative decay b —> sj serves as a probe of physics 
beyond the standard model (SM) . In the SM it proceeds 
at leading order through a loop diagram, making it sen- 
sitive to possible virtual contributions from as yet undis- 
covered particles. Because of parity violation in the weak 
interaction, the photon in fe — » 57 is predominantly left- 
handed, while it is right-handed in the charge-conjugate 
decay. The photon polarization can be determined in- 
directly through a measurement of time-dependent CP 
asymmetry in certain neutral decay channels. A non-zero 
asymmetry S due to interference between B^ mixing and 
decay diagrams is only present if both photon helicities 
contribute to both i3° and B^ decays fj . 5' is expected to 
be approximately —0.02 in the SM [l|,l4|, though hadronic 
corrections might permit it to be as large as ±0.1 [3]. Sev- 
eral new physics scenarios yield large values of the asym- 
metry; these include left-right symmetric models [l|, Q 
and supersymmetric models [5]. Because the SM asym- 
metry is small, any significant evidence of a large asym- 
metry would point to a source beyond the SM. 

We present an updated measurement of the time- 
dependent CP asymmetry in i?° -^ K^tt^^ based on the 
final BaBar data set of 467 x 10^ T{AS) -> BB decays 
collected at the PEP-II asymmetric-energy e+e" storage 
rings at SLAC. Previous measurements have been per- 
formed by BaBAR 6] and Belle 7] . Changes since BaBAR's 
last published result include doubling the data set, im- 
proved track reconstruction, better removal of back- 
ground photons from i:^ and 77 decays, better rejection 
of B'^ — + K*'^^ background ^], and an improved evalu- 
ation of the systematic uncertainties from non-signal B 
decays. At leading order in the SM, the CP asymme- 
tries of this mode do not depend on m(iir"7r'') [9|. How- 
ever, since the aforementioned hadronic corrections [3[ 
or new physics could introduce this dependence, we split 
the data into two parts: the K* region with 0.8 < 
m{K^TT^) < l.OGeV/c^ and the non-i^* region with 
1.1 < m{K^T:°) < 1.8GeV/c2. 



Time-dependent CP asymmetries are determined us- 
ing the difference of B^ meson proper decay times At — 
tsig — ttag , where f sig is the proper decay time of the sig- 
nal B° — > K'^tt^j candidate (Bgig) and itag is that of 
the other B (Btag), which is partially reconstructed and 
fiavor-tagged based on its daughter tracks. The At dis- 



tribution for Bsig decaying to a CP eigenstate is 



|At|/r 
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V±(At) = [1 ± S'sinfAmAt) T C'cos(AmAt)], 

4t 

(1) 
where the upper and lower signs correspond to Btag hav- 
ing fiavor B" and B° respectively, r is the B° lifetime, 
and Am is the B^-B^ mixing frequency. The C coeffi- 
cient corresponds to the direct CP asymmetry in decay, 
expected to be smaller than 1% in the SM |10 |. 

We evaluate our selection criteria with a detailed 
Monte Carlo (MC) simulation of the BaBAR detector fll| . 
using the EVTGEN generator "li\ and the GEANT4 pack- 
age [lj|. We require photon candidates to have energy 
greater than 30 MeV and the expected lateral shower 
shapes in the electromagnetic calorimeter (EMC). The 
primary photon from the B decay must be isolated by 
more than 25 cm from other charged and neutral clus- 
ters in the EMC. Primary-photon candidates that make 
a 7r° or ?7 candidate when combined with another photon 
in the event are discarded based on a likelihood formed 
from the diphoton mass and the energy of the second pho- 
ton. We select K^ — > tt+tt^ candidates from oppositely- 
charged tracks for which the probability of a geometrical 
vertex fit is greater than 0.1%, the tt+tt^ invariant mass 
is between 487 and 508 MeV/c^, and the reconstructed 
decay length is greater than 5 times its uncertainty. We 
select 7r° — > 77 candidates with invariant mass between 
115 and 155 MeV/c^ and energy greater than 590 MeV 
in the laboratory frame. For candidates in the K* re- 
gion we require jcos^i^. | < 0.9, where Ok* is the angle 
between the K^ and primary photon direction in the K* 
rest frame. 

To identify sign al decays we use the energy-substit uted 
mass mES = \/{s/2c^ + Po • VbY/EI - jpsp/c^ and 
the energy difference AB = E*^ — -\/s/2, where {Eq/c, po) 
and {Eb / c, ps) are the four-momenta of the initial e^e~ 
system and the B candidate, respectively, ^/s is the 
center-of-mass (CM) energy, and the asterisk denotes the 
CM frame. The distributions of signal events show a peak 
in these variables. We require 5.2 < toes < 5.3GeV/c^ 
and |AB| < 250 MeV. To reduce B+ -^ K*+-f back- 
ground, we reconstruct B+ -^ K*^ [KItt^)^ candidates 
subject to the same requirements as B° candidates, and 
veto events for which ?7ies(B+) > 5.27GeV/c^ and 0.8 < 
m{K2;Tr^) < l.OGeV/c^. To discriminate B decays from 
continuum e~^e~ -^ qq {q = u, d, s, c) background, we 
require |cos0|j| < 0.9, where 9g is the CM angle be- 
tween the B candidate and the e~ beam direction. We 
require the ratio of event-shape moments L2/B0 to be 
less than 0.55, where L, = J2j \Pj\\ cos 9* \\ p* is the CM 
momentum of each particle j not used to reconstruct the 
B candidate, and 6* is the CM angle between p* and the 
thrust axis of the reconstructed B candidate. After all 
selection criteria have been applied we find 10587 candi- 
date events, 16% of which have more than one signal B" 
candidate. In these cases we select the one with tt" mass 
closest to its nominal value UM , and if there is still an 



ambiguity, we select the one with the K'^ mass closest to 
its nominal value. We find an overall selection efficiency 
of 16%. 

For each reconstructed signal candidate we use the re- 
maining tracks in the event to determine the decay vertex 
position and flavor of -Stag • The latter is determined by 
a neural network based on kinematic and particle iden- 
tification information, the performance of which is eval- 
uated using a sample of fully-reconstructed, self-tagging 
hadronic 5° decays (iJflav sample) p^ . 

We determine the proper time difference between Bgig 
and Btag from the spatial separation between their decay 
vertices in the same way as our previous analysis and a 
similar BaBAR study of B° -^ K°Tr° ^. Because both 
the transverse flight length of the B° mesons and the 
transverse size of the interaction region are small com- 
pared to the B^ flight length along the boost direction, 
we are able to determine a decay vertex from the inter- 
section of the K'^ trajectory with the interaction region. 
We further improve the At resolution by 11% over what 
is obtained using information from the interaction region 
alone by refitting the T(4S') -^ B^B'^ system with the 
constraint that the average sum of decay times (tsig+^tag) 
be equal to twice the B lifetime with an uncertainty of 
\/2tb . Using MC simulation data we verify that this pro- 
cedure gives an unbiased estimate of Ai. We define events 
as having good At quality if each pion daughter of the K2; 
creates at least 2 hits in the silicon vertex tracker (SVT), 
and if the At uncertainty (jAt < 2.5 ps and | At| < 20 ps. 
About 70% of signal and background events pass these re- 
quirements. We split our data set and fitting procedure 
based on the At quality such that flavor-tagged events 
with poor At information do not contribute to the mea- 
surement of 5, but do contribute to the measurement of 
C, which can be determined solely through tagging. 

We extract signal yields and CP asymmetries using an 
unbinned maximum likelihood fit to tties, AE, L2/L0, 
tag fiavor, At, tJAt, and, in the K* region, m{K^'K^). 
Continuum and BB backgrounds are also modeled in 
the fit. We construct the likelihood function for each 
contribution as the product of one-dimensional probabil- 
ity density functions (PDFs). The signal PDFs in tties 
and IS.E are parametrized using the function 



f{x) = exp 



-{x-tif 



2cr2 + a{x - ^if 



(2) 



where /i is the mean, a the core width, and a a tail pa- 
rameter. The latter two parameters are allowed to be 
different on either side of the peak. The signal ra{K2:'K'^) 
shape is a relativistic Breit-Wigner, as nonresonant con- 
tributions in the K* region are negligible |17| . For con- 
tinuum rriES we use an ARGUS 18] function, while for 
l^E we use an exponential shape. The continuum and 
BB shape in ^^{K^tt^) is a Breit-Wigner on top of a lin- 
ear background. We parameterize the BB toes shape 
as the sum of an ARGUS function and a Gaussian with 
different widths below and above the peak, and the A_E 
shape as an exponential. The LiJLq shapes are binned 



PDFs, and in those the signal and BB components share 
the same parameters. All signal and BB PDF parame- 
ters are determined using simulated events, except for the 
flavor tag efficiencies, mistag probabilities, and At reso- 
lution function parameters, which are determined from 
the -Bfiav sample. The large number of continuum back- 
ground events in the fit determine the continuum PDF 
parameters. 

We obtain the At PDF for signal events and BB back- 
ground from Eq. [U accounting for the mistag probability 
and convolving with the At resolution function, which is 
the sum of three Gaussian distributions [l5| . The effec- 
tive GP asymmetries for the BB background, S ^^ and 

C„^, are fixed to zero in the fit, and we account for 
a possible deviation from zero in the systematic uncer- 
tainty. We verify in simulation that the parameters of 
the resolution function for signal events are compatible 
with those obtained from the Bflav sample. Therefore we 
use the Bfiav parameters for better precision. We fit the 
continuum MC At distribution and find that it is well- 
modeled by a prompt decay distribution consisting only 
of the At resolution function shape. The parameters of 
the continuum At PDF are determined in the fit to data. 
In the fit to the S" -^ K*^-^ candidate sample of 3884 
events we find 339 ± 24 (stat) signal events, Sk-^ — 
-0.03 ± 0.29 (stat) ± 0.03 (syst) and Ck*-i = -0.14 ± 
0.16 (stat) ± 0.03 (syst). We also find 19 ± 27 (stat) 
BB background events. In the range 1.1 < m{K2;'K'^) < 
1.8 GeV/c^ with 6703 events we measure 133 ± 20 (stat) 



signal events, Sj^o^o 



-0.78 ±0.59 (stat) ±0.09 (syst) 



and C^o^o^ = -0.36 ± 0.33 (stat) ± 0.04 (syst). We find 

167±49 (stat) BB background events in this sample. The 
linear correlation coefficient between Sk'j and Cx'-y is 
+0.050, while for Si^o^^o^ and C/^o^o^ it is +0.015. Fig- 
ure [1] shows signal-enhanced distributions for tties and 
AE created by cutting on the likelihood of the unplot- 
ted fit variables. 

We perform a cross-check because of the discrepancy 
between the projection of the fit model and the data in 
the non-K* region at low toes- A fit of the data sample 
with Toes > 5.22 GeV/c^ shows that the observed changes 
in S and C are consistent with statistical fluctuations, so 
the signal is not significantly affected. Additionally, we 
verified that the slope of the to-es background shape is 
not correlated with the other fit variables. MC simu- 
lations of common BB backgrounds, including the final 
states K'^ttttj and K'^mnr'^ , do not show any rising struc- 
ture at low m-Es- 

Figure [5] shows the background-subtracted distribu- 
tions of At in the K* region, obtained with the sPlot 
event weighting technique Jl!| . We show an sPlot of the 
m(i^g7r'^) spectrum in Fig. [3] 

Using an ensemble of simulated experiments generated 
from the fitted likelihood function we find no bias in the 
K* region and a spread in 5" and C consistent with the 
statistical uncertainties. In the non-_ftr* region we find a 
bias of -0.06+0.03 on S and a spread in C larger than the 
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FIG. 1: Signal-enhanced distributions for rriES (top) and 
AB (bottom) for the K* region (left) and the non-A"* region 
(right). We show the fit result (solid line) and PDFs for signal 
(long dashed), continuum (short dashed), and BB (dotted). 
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FIG. 2: sPlot (see text) of At in the K* region (left) and the 
non-Zc" region (right), with Stag tagged as B (top) or B 
(center), and the asymmetry (bottom). The curves are the 
signal PDFs. 



statistical uncertainty reported by the fit. These effects 
are due to a measurement that is close to the physical 
boundary of S"^ + C^ < 1 , and they disappear if we gen- 

= Cxo-^o^ = 0. The bias 



erate the ensemble with Sj^o^o^ 



on 5 



KgTrS 



is evaluated with several ensembles of simu- 
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FIG. 3: sPlot (see text) of m(Kl-!T" 



the statistical uncertainty on C we take the ensemble's 
root-mean-square width of 0.33 instead of the 0.29 un- 
certainty determined by the fit to data. 

Systematic uncertainties associated with our knowl- 
edge of the beam spot position and possible SVT mis- 
alignment are determined by varying the beam spot and 
SVT alignment parameters in MC. We bound the effects 
of uncertainties in the At resolution function due to the 
vertexing method with a study from BaBAR's i?" -^ ^3""° 
analysis [20]. Resolution function differences between 
data and MC in control samples of B'^ -^ J/ipK^; decays, 
in which the J/ip vertex information is ignored, lead to 
differences in S and C that we take as systematic uncer- 
tainties. Uncertainties from doubly-Cabibbo-suppressed 
(DCS) decays of the Stag are included as in Ref. [l5| . 

We evaluate uncertainties due to the vertex reconstruc- 
tion procedure and possible correlations among the ob- 
servables with an ensemble of simulated experiments cre- 
ated by generating background events from the PDFs 
and embedding signal events from the full MC simula- 
tion. No significant bias is observed in the K* region, 
and we bound uncertainties by the precision with which 
the potential bias is measured. In the non-_R'* region, no 
bias is observed in the signal MC. 

Uncertainties due to limited knowledge of the fixed pa- 
rameters in the fit are evaluated by varying them within 
their uncertainties. We evaluate differences between data 
and MC in the signal shape by fixing the background pa- 
rameters to those determined in the fit to data and float- 
ing the signal parameters separately for each observable. 



We evaluate the effect of S" - and C - by varying 

BB BB •' JO 

them over a range determined by the composition of the 
BB background samples and CP asymmetry measure- 
ments in the PDG listings. The systematic uncertainties 
are summarized in Table [H 

In summary, we have measured the time-dependent CP 
asymmetry in 5° -^ K'^ir'^j decays using the full BaBAR 
data set recorded at the T{4S) resonance. We find 



lated events in which the generated 5^0 ^,.0-, is varied. For 



Ck'j = 



-0.03 ± 0.29 (stat) ± 0.03 (syst), 
-0.14 ± 0.16 (stat) ± 0.03 (syst), 
-0.78 ± 0.59 (stat) ± 0.09 (syst), 



TABLE I: Summary of systematic uncertainties. 





K* 


Region 


non-K* 


Region 


Source 


AS 


AC 


AS 


AC 


Beamspot 


0.007 


0.002 


0.007 


0.002 


SVT Alignment 


0.010 


0.010 


0.010 


0.010 


Resolution Function 


0.011 


0.018 


0.011 


0.018 


Bias Uncertainty 


0.015 


0.009 


0.028 


0.016 


PDF Uncertainty 


0.015 


0.013 


0.060 


0.019 


S-'^'^i and C'"'i 

BB BB 

DCS Btag Decays 


0.008 
0.001 


0.002 
0.015 


0.060 
0.001 


0.018 
0.015 


Total 


0.028 


0.030 


0.091 


0.040 



Ciio^o^ = -0.36 ± 0.33 (stat)± 0.04 (syst). 
The measurement in each 'm{K'^Tr'^) region is consistent 



within uncertainties with the predictions of the standard 
modcL 
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